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The Orientation of the Reconnection X-line 
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We propose a criterion for identifying the orientation of 
the X-line when two regions of plasma with arbitrary den- 
sities, temperatures, and magnetic fields undergo reconnec- 
tion. The X-line points in the direction that maximizes the 
(suitably-defined) Alfven speed characterizing the reconnec- 
tion outflow. For many situations a good approximation is 
that the X-line bisects the angle formed by the magnetic 
fields. 



1. Introduction 

Reconnection is the dominant process by which energy 
is transferred from the magnetic field to the thermal and 
bulk motions of the particles in collisionless plasmas such 
as the magnetosphere and the solar corona. Both theoret- 
ical models and numerical simulations of reconnection usu- 
ally consider highly symmetric cases, e.g., the merging of 
two plasmas that are identical except for their anti-parallel 
fields, where symmetry considerations dictate the reconnec- 
tion plane and the orientation of the X-line (the normal to 
that plane). Realistic configurations are often more com- 
plex, as for instance at the magnetopause where a low- 
density, strong-field plasma (the magnetosphere) merges at 
an arbitrary angle with a high-density, weak-field plasma 
(the magnetosheath). 

Sonnerup [1974] argued that in such complex systems the 
orientation of the X-line is fixed by requiring that currents in 
the reconnection plane vanish, and hence, by Ampere's Law, 
that the guide field (the magnetic component parallel to the 
X-line) in the two plasmas be equal. However this choice has 
the peculiar consequence that there are some magnetic field 
configurations for which reconnection cannot occur because 
the reconnecting components of the field have the same sign. 
A further concern arises from the observation that when a 
thermal pressure gradient exists at an X-line the guide field 
must have spatial variations if the system is to be in total 
pressure balance. Since there is no a priori reason for assum- 
ing thermal pressure gradients vanish at X-lines this calls the 
primary motivation for Sonnerup's choice into question. 

We propose a different criterion: reconnection occurs in 
the plane in which the outflow speed from the X-line (given 
by an appropriately-defined Alfven speed) is maximized. 
With this choice reconnection can occur between any plas- 
mas in which the magnetic fields are not exactly parallel. 
Reassuringly, the orientation of the X-line also reduces to 
the expected result in symmetric cases. 

2. Definition of Coordinates 

It is particularly important for this problem to define the 
coordinates carefully. Consider two regions of plasma each 
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with number density rij , temperature Tj , and magnetic field 
Bj, where j — 1,2. Assume that the two regions are sepa- 
rated by a planar discontinuity through which no magnetic 
field passes and define a coordinate system in which the x 
and z axes lie in the discontinuity plane and the y axis is 
perpendicular to it. 

Without any further constraint the X-line could, in prin- 
ciple, point in any direction in the x — z plane. Each different 
X-line orientation implies a different reconnection plane with 
different components of the field reconnecting and a different 
reconnection rate. We want to find the X-line orientation 
for which reconnection is fastest. To do so it is most con- 
venient not to consider a fixed coordinate system in space 
but rather to define our coordinates with respect to the di- 
rection of the reconnection X-line (the z axis) and the plane 
of reconnection (the x — y plane) and to rotate the fields 
about the y axis. This rotation intermixes the guide, z, and 
reconnecting, x, components of the fields and changes the 
reconnection rate. The GSM equivalents of our coordinates 
at the magnetopause are (x,y,z) —* (z,x,y)csM- 

Without loss of generality we specify the orientations of 
the fields by defining 9 to be the angle between the fields on 
either side of the discontinuity (also called the shear angle) 
and a as the angle Bi makes with the z axis (see Figure 1). 
To make the problem well-defined we limit the ranges of the 
angles: < 9 < it and < a < 9. 

The unknown parameter is a and varying a at fixed 9 
changes the relative orientations of the fields with respect 
to the X-line. According to Sonnerup's argument a is the 
solution of the equation Bi cos a = B2 cos(# — a). We claim 
that the proper choice is instead the a that maximizes the 
outflow speed and the rate of reconnection. 

As an example, consider a system with 9 — n, Bi = 
— B2, ni = ri2, and Ti = Ti. These parameters describe 
anti-parallel reconnection and symmetry suggests that a — 
tt/2 — 9/2. Adding a constant guide field will change 9 but 
should keep a = 9/2. Sonnerup's criterion gives the ex- 
pected results in these cases and ours, as will be seen, does 
as well. For other parameters, however, the two differ. 

3. Determining a 

The rate at which magnetic field lines reconnect directly 
varies with the speed at which they flow toward the X-line. 
Continuity suggests that the speed of this inflow is propor- 
tional to the speed of the field lines' outflow, with a constant 
of proportionality that depends on the detailed physics of 
the reconnection (e.g., the aspect ratio of the diffusion re- 
gion). For our purposes the details of the dependence do 
not matter; the crucial point is that as the outflow speed 
increases the reconnection rate does as well. 

Since the outflow is driven by the motion of magnetic field 
lines it must be related to some Alfven speed; for symmetric 
anti-parallel reconnection it is the speed calculated from the 
asymptotic field and density. Defining the appropriate out- 
flow speed in the general case is more complicated. We find 
that it depends on the fields and densities in both plasmas 
as well as the angles 9 and a. Hence, the inflow speed and 
reconnection rate depend on these quantities as well. 

3.1. Constructing the outflow speed 

Consider the situation shown in Figure 2. The left side of 
the box is the X-line where the plasma velocity and in-plane 
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magnetic field are assumed to vanish. The plasma above 
the current sheet has mass density pi — mini, where mi is 
the average mass, and in-plane magnetic field B Xt i; below 
the current sheet these values are pi and fix, 2. The plasmas 
flow into the current sheet with speeds vi and «2- Within 
the current sheet they mix in some proportion, resulting in 
a plasma of density po, and accelerate downstream, dragged 
by the tension of the reconnected magnetic field. At the 
right-hand edge of the sheet the plasma reaches its outflow 
speed vo and the in-plane field has a magnitude B V] o- We 
assume the system is in a steady-state and proceed to cal- 
culate Ho. 

Applying conservation of mass to the box gives 



pov Ay = piviAx + P2V2AX . 



(1) 




The out-of-plane electric field E z is, according to Fara- 
day's Law, spatially constant in a 2-D steady-state sys- 
tem and, asymptotically, must be given by the MHD re- 
sult E = —v X B/c. Equating the values at the inflow and 
outflow edges of the current layer gives 

ViB x< i = v 2 B Xi2 = v Byfi ; (2) 

combining equations (1) and (2) yields an expression for po 
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Within the current layer the dominant terms in the x compo- 
nent of the momentum equation are advection and magnetic 

tension: a 1 an 

dv x 1 dB x 

PoV *-te = ^ B »-dy-- (4) 

We assume B x varies piecewise-linearly across the current 

layer and rewrite this equation as 



d 2 1 By B x 1 + B x 2 

— v x = — — — 

ox 47r po Ay 



(5) 



After integrating with respect to x along the current layer 
we have 

2 _ 1 B y , Ax 
V °-^—Ay- {Bx - 1+Bx - 2h (6) 
where J B y dx = B y fiAx. Combining equations (3) and (6) 
gives the outflow speed 
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Equation (7) exhibits the necessary symmetry between the 
two sides, reduces to the usual result, Vq = B^/A-Kp, when 
pi = P2 and Bi x = B2 X , and goes to zero, as expected, 
when either density is large or either field vanishes. This re- 
sult was independently derived in a slightly different context 
by Cassak and Shay [2007]. 

In terms of the angles defined in Figure 1 the outflow 
speed is 



Figure 1. Definition of the coordinate system. The plas- 
mas meet in the plane shown, with one extending out of 
the page and the other into the page. The shear angle 
9 and the directions of the x and z axes are fixed; a is 
unknown. 
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Figure 2. Cartoon of the current layer during asymmet- 
ric reconnection. The X-line is at the left of the box and 
the asymptotic outflow is at the right. All speeds and 
magnetic field strengths are assumed to be positive. 
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Figure 3. Plot of a versus the shear angle 8 for the case 
b — 2, r — 0.5. Note that for Sonnerup's solution there 
is no reconnection when 9 < 60°. 
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According to our previous argument the condition dv 2 . /da = 
defines the orientation of the X-line. 

3.2. Maximal Value 

Although the operations required to find an expression 
for a are straightforward, the actual calculations are a bit 
tedious. Before presenting the result, we make some obser- 
vations 

1. v 2 (a = 0) = v 2 (a = 9) = 0. Since v 2 > the im- 
plication is that Vq has at least one maximum in the range 
< a < 6. We strongly suspect, but have not been able to 
prove, that there is only one maximum. 

2. B\, £?2, pi, and p2 are independent variables but will 
only enter the result through the two dimensionless ratios 
b — B2/B1 and r = pi/ p\- Hence a is a function of only 
three parameters: 6, b, and r. 

The solution for a is the root of the equation 

= r sin 2 a[sin(0 - 2a) - 6sin(20 - 2a)] 

+ b sin 2 (6 - a) [sin 2a + b sin(0 - 2a)] (9) 

subject to the constraint < a < 6. By defining (f> — 6/2— a, 
4> + — 6/2 + <f>, and 0_ = 6/2 — <f> equation (9) can be written 
in the symmetric form 

= r sin 2 0_ [sin 2<p — b sin 2(j> + ] 

+ 6sin 2 0+[fosin20 + sin20_]. (10) 

Although equation (9) must, in general, be numerically 
solved for a, exact solutions are possible in some special 
cases 

1. 6 = 7r (anti-parallel reconnection) . In this case a = 
6/2 = 7r/2, independent of the values of b and r. 

2. r = 1 (pi — p2)- Regardless of b the maximal value 
occurs for a — 6/2. 

3. b > |1 - l/r\ or 6 « 1 1/(1 - r)\. Again the re- 
sult is a — 6/2. The two limits are complementary in the 
sense that the system is symmetric under the substitutions 
b — > 1 /&* , r — > 1 /r« , a — > 6* — a* . 

The last example suggests that a = 6/2 is a good approx- 
imation to the exact solution of equation (9) whenever the 
density ratio is not too much different from 1. Numerical 
trials bear this out, as can be seen in Figure 3 which shows 
results for b — 2, r — 0.5. 

Interestingly, since the outflow speed does not directly 
depend on the temperatures or average masses of the plas- 
mas, neither does a (or, equivalently, the X-line orientation) . 
There is an indirect constraint, however, because the system 
must also be in total pressure balance, 



B 2 B 2 
niTi + -i = n 2 T 2 + . 
87r 8n 



(11) 



if our assumption of steady-state reconnection is to be valid. 
If the temperature and the average mass are equal in the re- 
connecting plasmas then equation (11) relates b and r to the 
plasma 



4. Discussion 

Establishing the system's orientation is an important part 
of the interpretation of spacecraft observations. Beginning 
with the basic magnetic field data the well-known technique 
of minimum variance analysis determines the direction nor- 
mal to the current sheet (the y axis in our coordinates). 
Determining the direction of the X-line, either through Son- 
nerup's criterion (see, for example, Phan et al. [2006]) or 
through equation (9), fixes the geometry of the reconnection, 
provided only that the system has weak variations along the 
direction of the X-line. This information is particularly im- 
portant for those measurements that are to be compared 
to theoretical models and simulations of reconnection, as 
will be the case for the upcoming Magnetospheric Multi- 
scale Mission. 

Our proposed criterion can be checked with numerical 
simulations. Since our argument does not depend on the 
detailed physics of the reconnection, only that the recon- 
nection rate varies with the Alfven speed, even MHD codes 
that do not correctly describe fast reconnection should suf- 
fice. However such simulations must take care not to impose 
a reconnection plane a priori by, for example, not being fully 
three-dimensional. 

We emphasize that although we have attempted to calcu- 
late the direction of the dominant reconnection X-line in a 
general current layer in this paper, there are several possibly 
important effects that have been neglected. First, we cannot 
exclude the possibility that reconnection may proceed simul- 
taneously at different surfaces and that, as a consequence, 
the current layer might become fully turbulent [Galeev et al, 
1986]. Second, effects that preferentially suppress reconnec- 
tion for some X-line orientations are a possible complication 
that we have ignored. Swisdak et al. [2003] showed that a 
thermal pressure gradient across the current layer drives dia- 
magnetic drifts that convect the X-line. As the drift speed 
approaches the Alfven speed the reconnection can be com- 
pletely suppressed. Since the magnitude of the drift varies 
with the angle a, the X-line orientation in such systems 
may be determined by a trade-off between maximizing the 
outflow Alfven speed and minimizing the diamagnetic drift. 
Other effects, e.g., shear flows in the reconnecting plasmas, 
could have similar consequences. 

Finally, equation (9) determines the local orientation of 
the X-line based on the parameters of the reconnecting plas- 
mas. But what happens at, for instance, the magnetopause 
where the shear angle can vary with location due to the 
combined effects of the dipole tilt of the terrestrial field, the 
direction of the interplanetary magnetic field, and the cur- 
vature of the interface? Both the orientation of the X-line 
and the reconnection rate will then vary with location with 
unknown effects on the global configuration of the reconnec- 
tion. One possibility is that local maxima in the reconnec- 
tion rate will seed vigorously growing X-lines that propagate 
outwards [Huba and Rudakov, 2002], perhaps occasionally 
shifting directions to merge with other reconnecting regions. 
Depending on the external conditions and length of time the 
system remains in a steady-state it may have either a few 
or many simultaneously reconnecting X-lines. 
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If desired the condition of equation (9) can be re-written in 
terms of fii and 02- 
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